Controllable half-vortex lattices in an incoherently pumped polariton condensate 
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We show how the transition between synchronized and desynchronized states of a spinor polariton 
condensate can be used to drive a transition between stationary vortex lattices and half-vortex 
lattices. This provides a way to control polariton spin textures by a combination of pump spot 
profile and applied magnetic fields. To do this, we extend the model of non-equilibrium spinor 
condensates to include relaxation, and study how this affects the desynchronization transition. We 
discuss how the pattern formation can be explained by behavior of the homogeneous system. 
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Among the various topological defects that may exist 
in superfluids and superconductors, quantized vortices 
receive particular theoretical and experimental attention. 
These localized nonlinear structures are nontrivial low 
energy excitations of the system, and so knowledge of 
the core structure, topology, symmetry, dynamics and 
interactions of vortices reveals details of the low energy 
properties of the system. As well as conventional sin- 
gular quantized vortices seen in ^He-II superfluids, one 
component Bose-Einstein condensates, and s-wave super- 
conductors, there also exist more exotic quantized vor- 
tices, such as non-singular vortices, continuous vortices 
and composite vortices Non-singular vortices have a 
hard core set by the coherence length but with a non-zero 
order parameter within the core. Such order parameter 
textures have been explored in various systems, includ- 
ing multi-component cold atomic gases [2| , nonlinear op- 
tics 0], superfluid ^He-A[H and Sr2Ru04[l. However, 
the control of these structures and, in particular, their 
nucleation from conventional singular vortices has been 
elusive. In this letter we propose how these structures 
can be controUably created and manipulated in polari- 
ton condensates by using magnetic fields. We make use 
of the non-equilibrium properties of the polariton con- 
densate, and so illustrate how a nonequilibrium conden- 
sate may show behavior not present in equilibrium spinor 
condensates like cold atoms or ^He-A. 

Microcavity polaritons are quasiparticles that result 
from the hybridization of excitons and photons confined 
inside semiconductor microcavities. At low tempera- 
tures, they can condense, and polariton condensation has 
been the subject of much recent work (see Q for a re- 
cent review). Polariton condensates differ from conden- 
sates of ultracold atomic gases or helium in that polariton 
condensates can be non-equilibrium systems: the photon 
component decays, and so new polaritons must be in- 
jected into the system to maintain a steady state. Hence, 
polariton condensates can be thought of as sustained 
nonequilibrium systems with nontrivial spatiotemporal 
properties. Spatiotemporal pattern formation has been 
extensively studied in systems such as fluid and plasma 



dynamics, reaction diffusion systems, biological struc- 
tures, and neural networks. Order parameter equations 
describe the structure and dynamics of pattern formation 
in these systems, and the order parameter equation for 
the polariton system (and closely related nonlinear opti- 
cal systems [61) is the complex Ginzburg-Landau equa- 
tion (cGLE) [7|. Polariton condensates are well suited 
for studying pattern formation since the source term can 
be made to vary in space (and time) and thereby con- 
trol pattern formation. Recent experiments on polari- 
tons have begun to explore this possibility, by studying 
phenomena such as quantized vortices 0] , soliton propa- 
gation , and patterns in onc-dimcnsional samples flo| . 

As polaritons have a polarization degree of freedom, 
they can support unconventional quantized vortices, such 
as non-singular vortices where only one polarization com- 
ponent has vortices [111. Such polariton "half- vortices" 
have been seen in experiments 12], but there are ques- 
tions about the conditions required for them to remain 
stable and independent (l3|. Experimental investigation 
of their stability is hard because observations of polari- 
ton vortices have generally relied on disorder to pin vor- 
tices and half- vortices could be stabilized by disor- 
der induced strain splitting of polarization states. The- 
oretical alternatives to these disorder-localized vortices 
have been proposed, such as rotating lattices in harmon- 
ically trapped incoherently pumped condensates • Di- 
rect observation of such rotating vortex lattices is how- 
ever hard with current approaches to observing vor- 
tices 15,1^. There are also proposals to create half vor- 
tices using parametric 17, 13] or coherent pumping [l9| . 
None of these schemes combine the necessary features to 
experimentally investigate the stability of half vortices in 
a clean incoherently pumped system. This is the task we 
address in this letter, by showing how interplay between 
the synchronization-desynchronization transition of the 
homogeneous system ^Q, [2lj and interference between 
polaritons originating from separate incoherent pumping 
spots can produce non-rotating lattices of vortices, where 
a transition between separate and locked half- vortices can 
be controlled by an applied magnetic field. 
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We start by introducing the cGLE for a spinor polari- 
ton condensate. This includes a term describing pump- 
ing, and for stabihty this must be nonhnear, so that it 
saturates as the condensate density increases. In a sin- 
gle component condensate, the simplest such nonlinear- 
ity is to write, ihdtip ~ 1(7 — r|?/;p)'0. However, as dis- 
cussed recently [l^], adding an energy dependent gain 
term ihdtip ~ «(7 — Tj-^p — ir]hdt)'ip can be important 
to correctly describe pattern formation where energy re- 
laxation plays an important role. The (dimensionless) 
parameter rj is the rate at which gain decreases with in- 
creasing energy. In the spinor case, a further possibility 
emerges: there can be a cross-polarization term whereby 
the gain of one component may be affected by the density 
of the other component. These various nonlinearities can 
all be understood in terms of transfer between an exci- 
tonic reservoir and the condensate, driven by the lack of 
chemical equilibrium (23j . Nonlinearity arises both from 
depletion of the reservoir, and from shifts of the conden- 
sate chemical potential. We have therefore included all 
three types of possible nonlinearity, and studied how they 
affect the behavior. The spinor cGLE is thus: 
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+iP(r,?A±)-iK]V'± +f}||V'=F, (1) 



where ± stand for left- and right- polarizations. We first 
discuss the parameters describing energetics, and then 
turn to those related to the pump configuration and loss. 

Since the nonlinear terms can be rescaled by a change 
of density, the only important parameter for the interac- 
tions is the ratio Ui/Uq which we take as Ui/Uq = 0.55 
(see e.g. [Ilj). This implies a weak attraction between 
opposite polarizations of polaritons. The parameter n± 
represents a magnetic field perpendicular to the 2D po- 
lariton condensate, while r2|| represents either an in- 
plane field, or the intrinsic splitting of linear polariza- 
tions present in most microcavity samples. We consider 
values based on the recent experiment in GaAs[24|, for 
which a range < ^± < 50/ieV was achieved, result- 
ing in a transition from linear to circular polarization. 
The splitting riy was only constrained to be less than the 
linewidth. We choose a value fiy = 3/ieV consistent with 
this, which, as discussed below, reproduces results for the 
homogeneous system that may be compared to experi- 
ment. Finally, the polariton mass only controls the over- 
all length scale io = ^ h? /2mK\ where necessary we give 
lengths in terms of this scale which would be €0 — 1.95/xm 



for m = 10 



taking decay rate k = O.lmeV. 



We propose to form vortex lattices by pumping in three 
localized spots equidistant from each other. This means 
that the pump term P{r,'ip±) is taken to be zero out- 
side the pump spots. Inside each pump spot, includ- 
ing the nonlinearities discussed above, it has the form 



pumping at 2.5 times threshold, i.e. 7 = 2.5k. The 
parameter rj is not well constrained, so we consider a 
range of values. For the other nonlinearities, we take 
Ts/Un — 0.3,Tx/Uo = 0.1, based on values discussed 
in j2l| . The term Fx, describing a cross-spin nonlinear- 
ity is new, and we do not have an independent estimate 
of its value. However, as long as F^; < F^, we find that 
the effect of increasing F^^/Fj is similar to the effect of 
increasing 77. If F^, = F^, the pumping process becomes 
spin-blind, and so the polarization state is fixed purely 
by the energetics. 
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FIG. 1. Vortex lattices due to polariton flow from three pump- 
ing spots. Panel (a) shows a 40^o x 40^o region, illustrating 
the location of the pumping spots and a smaller 2£o x 2£o re- 
gion marked by a (red) square. Panel (b) shows the stream- 
lines and density profile of the majority component at il± = 
2.4^eV in the synchronized regime. (The minority streamlines 
are locked and thus equivalent). Panel (c) Show the stream- 
lines for both polarizations (Red (gray)- majority, Black - 
minority) in the desynchronized regimes for i}± = lOO/xeV 
All figures are for 7 = 2.5k, — 3fieV,ri — 0.1. 

If the pumping spots are irregularly spaced, then the 
generated currents interfere and form vortex pairs which 
may lead to turbulent state of interacting vortices as 
was discussed in [2^. However if the spots are placed 
sufficiently symmetrically then stationary vortex lattices 
form. We performed series of numerical experiments 
pumping at three spots of radius R = 5£o centered at 
(0, 20)4, (±10^3, 10)£o (see Fig. [IJ. The numerical in- 
tegration of Eq. ^ was done by the fourth order finite 
differences in space and the fourth order Runge-Kutta 
in time. Vortex lattices form as long as the flux leaving 
one spot is enough to reach the adjacent spots, and thus 
cause synchronization; this is very similar to what was 
found for polariton condensates trapped in disorder po- 
tentials [20|, however here we have no trapping potential. 
Vortex lattices are thus visible only for large enough j/k 
and small enough 77. Fig. [T] is for parameters where syn- 
chronization between spots just occurs, and as a result, 
the lattice is not perfectly hexagonal; for larger a, lat- 
tices are more hexagonal (not shown). Synchronization 
of pump spots, and the consequent vortex lattices, can 
also survive some degree of imperfection in size and loca- 
tion, 10% variations can be tolerated for the parameters 
used in Fig. [1] 

For conditions where a vortex lattice occurs, we then 
vary the detuning n± and observe the evolution of vortex 
lattices. As 1 increases, there is a transition from an 
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elliptically polarized condensate to a state where the two 
circular polarization components desynchronize, produc- 
ing two separate condensates, with two separate chemical 
potentials (frequencies). In the synchronized regime, in- 
teger vortex lattices are formed with vortices of the same 
circulation of the two polarizations sitting on the top of 
each other. In the desynchronized regime, the lattice of 
the majority component survives, but a number of possi- 
ble behaviors for the minority component can occur, de- 
pending on the values of 7/K, ij, ft±. In Fig.fljc), the mi- 
nority component has a vortex-free state. There can also 
be cases where a vortex lattice survives for the minority 
component, but because the polarization components are 
not synchronized, the position of these minority vortices 
can oscillate in time. For large j/k, these oscillations are 
small, and only a slight shift of the minority lattice is vis- 
ible. However, for 7 = 2.5k as shown in Fig.[ljc), there is 
a clear transition to a half-vortex lattice. At much larger 
, a transition back to a synchronized state is seen, but 
for 77 = 0.1, this is far beyond the range of ft± achievable 
in experiment. 

Having illustrated how the desynchronization transi- 
tion has dramatic effects on the vortex lattices, we derive 
analytic approximations for the period of the vortex lat- 
tice, and the critical field for the desynchronization tran- 
sition, by using a generalized local-density approxima- 
tion (considering also the phase gradients), and making 
use of results for a homogeneously pumped system. The 
period of the vortex lattice can be understood by finding 
the wavelength of the phase variation at points far from 
a pump spot. If the pump spot is large then currents 
can be neglected within the pump spot, and the homoge- 
neous system used to predict whether the system is syn- 
chronized, and the values of chemical potential. When 
synchronized, idtip± = fjL'>p±, and so Eq. ([1]) becomes a 
pair of time-independent differential equations. To solve 
these, it is convenient to write ip± = y/])±e^^'^^^ In 
the synchronized case 6 is constant so the lattice spac- 
ing is set by the velocity u = hV(j)/m. Far away from 
the pump spot, the pump term vanishes, and the density 
decays and so interaction terms U p± can be neglected. 
With these approximations the real and imaginary parts 
of Eq. (IT} become: ^~^m\\i\'^^^l±^ — ^ p±/ cos{6) 
and |fi.V • (/5±u) = —Kp± ± 0|| ^ p^p^ sin(6'). These 
equations are solved by sin(0) = 0, p+ cx p- and thus 
^m|up — p- ^ (f^l + When desynchronized, a 

general solution is more complicated, as Q and p-\-/ p- are 
time dependent, and one must find the average (|u±p) 
in terms of {p±) ■ However, if p- ^ p^ in the desyn- 
chronized region, then the real part of Eq. ([T]) for the ma- 
jority component becomes ^m(|u_p) = (^_) -1-57^. The 
hexagonal vortex lattice resulting from the interference of 
three such currents has a spacing I = (2/3v^)(/i/m|u|). 
For the parameters of Fig. [1] this gives / = 1.4£o, con- 
sistent with the observed lattice. For larger rix,51||, the 



lattice spacing is found to vary considerably with n±_ ; we 
have checked that the form discussed here matches the 
numerical simulations well. 
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FIG. 2. Panel (a): boundary of desynchronized region as a 
function of field Q± and relaxation rate rj for homogeneous 
pumping. Solid lines are exact numerical results; gray dashed 
lines are the approximate asymptotes given in the text. Panel 
(b): polarization degree vs. field Six for the homogeneous sys- 
tem at 77 = 2. Lines indicate the steady state polarization for 
the synchronized regimes. The points and error bars indicate 
the time-average, and envelope of polarization oscillations. 
(NB the envelopes are normalized by the average, rather than 
instantaneous, intensity.) 

We next discuss the critical fields required to switch 
between synchronized and desynchronized lattices, as 
shown in Fig. [2ja). The line in Fig. [2] is calculated for 
homogeneous pumping; numerical simulations with in- 
homogeneous pumping also respect this boundary, i.e. 
the critical r2j_ for desynchronization is rather robust to 
inhomogeneity. Analytical estimates of the critical fields 
ri_L, help to reveal how they depend on the system pa- 
rameters. The first transition on increasing fl± is from 
the elliptically polarized (synchronized) phase existing 
for small il± to a desynchronized state. This transition 
was discussed in Refs. [23,|2l|. Ifrij_,ri|| < 2Ui{p++p-) 
(which is the case for the current parameters), then 
one may approximate the coupled equations for 9, p± 
(with the parametrization tlj± = ^p^e*'"^^^/^^ intro- 
duced above) by using standard results for Josephson 
oscillations in the Josephson regime [2^. The transition 
thus occurs at n± = ilci — 2Uiil^^/{Ts — T^)- 

If 77 7^ there is another transition at = flc2 back 
to a synchronized state. The critical field flc2 can be 
estimated analytically by linear stability analysis of the 
synchronized state at 0_l > Qc2- Since p_ 3> p+ in this 
large fl± limit, one may show that fiuctuations of — 
p_ relax quickly. The remaining part of the linearized 
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Eq. (m can be written in terms of = ^is: 

(1 + rf)hdtC = (7cfr - i^esK - ^\\P-iv + (2) 

where 7eff = 2(7 - k) - p„[r, + + 2'ij{Uo - Ui)] and 
j/off = 2Q± + [77(rs — r^) — 2Ui]p-. Instability of the syn- 
chronized state occurs when ^^s becomes negative, thus 
to find ilc2 one thus needs at large rt±. Since 

p_ ^ p^, one may simplify Eq. ([T]) for ip^ considerably 
to give p_ ~ (7 — K + rjH.^)/ (rjUo + Tg) hence one has 



(7- 



77 Ts+T, + 2r,{Uo-Ui)' 



(3) 



As the relaxation rate 77 increases flc2 becomes smaller, 
and eventually the desynchronized regime will vanish. 

Let us now briefly consider the relation of the exper- 
imental results of Larionov et al. [2J] to the polariza- 
tion desynchronization transition. In their experiment, 
a transition was seen such that for small a single 
frequency elliptically polarized condensate was observed, 
and above a critical value of fij, , a transition to a Zeeman 
split condensate, with separate frequencies for left- and 
right-circular components was observed. While it was 
noted in [gj] that some features of the experimental re- 
sults clearly required aspects of non-equilibrium physics, 
the transition was interpreted as a signature of the equi- 



librium transition to a circular condensate [27[. However, 
a fragmented condensate (i.e. having multiple conden- 
sate frequencies) is not expected in equilibrium, but is 
precisely the signature expected for the desynchroniza- 
tion transition. While this identification would require 
further experimental corroboration, we have chosen our 
value of r2|| = 3//eV so the desynchronization transition 
matches this experimental transition. Figure ^h) illus- 
trates the behavior in the desynchronized regime, show- 
ing the amplitude of polarization fluctuations. The time 
averaged polarization degree shown in that figure evolves 
smoothly with field, and so provides a poor distinction 
between synchronized and desynchronized behavior. One 
may however note that both f2c2, and the equilibrium 
critical field j27j increase with density. In contrast, fid, 
where splitting of the condensate frequencies should oc- 
cur due to desynchronization is independent of density 
(for small enough 77). This may provide a clear way to 
distinguish the nature of the transition observed experi- 
mentally 0. 

In conclusion, we have explored the desynchronization 
transition in a non-equilibrium spinor polariton conden- 
sate, and shown how this transition can be used to drive 
a system between lattices of integer and half vortices. 
Using the parameters of recent experiments 2J] , we show 
that the conditions for the desynchronization transition 
seem likely to be achievable, and note that aspects of 
the experimentally observed behavior may even suggest 
that the transition already seen is that to a desynchro- 
nized state, rather than to a circularly polarized state. 



When the ability to control polarization via magnetic 
field is combined with the possibility to form static vor- 
tex lattices via symmetric pumping spots, our proposed 
scheme illustrates how one may engineer and control ex- 
otic lattice spin textures deterministically using incoher- 
ently pumped polariton condensates. 
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